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The  axial  compressive  mechanical  response  of  substrate-supported  carbon  nanotube  (CNT) 
arrays  with  heights  from  35  to  1200  nm  is  evaluated  using  flat  punch  nanoindentation  with 
indentation  depths  to  200  fim.  The  compressive  behavior  is  consistent  with  that  of  an  open¬ 
cell  foam  material  with  array  height  playing  a  role  similar  to  that  of  occupation  density  for 
traditional  foam.  Mechanical  yielding  of  all  arrays  is  initiated  between  0.03  and  0.12  strain 
and  arises  from  localized  coordinated  plastic  buckling.  For  intermediate  CNT  array  heights 
between  190  and  650  pm,  buckle  formation  is  highly  periodic,  with  characteristic  wave¬ 
lengths  between  3  and  6  pm.  Buckle  formation  produced  substantial  force  oscillations  in 
both  the  compressive  and  lateral  directions.  The  compressive  elastic  modulus  of  the  arrays 
is  obtained  as  a  continuous  function  of  penetration  depth  and  attains  a  value  between  10 
and  20  MPa  for  all  arrays  during  mechanical  yield.  A  qualitative  model  based  upon  concepts 
of  cellular  foam  geometry  is  advanced  to  explain  the  observed  CNT  buckling  behavior. 
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1.  Introduction 

Aligned  CNT  arrays  have  received  significant  attention  for 
potential  use  in  applications  such  as  thermal  [1-7]  and 
electrical  [8-11]  interface  materials,  electrical  interconnects 
[12-14],  and  composite  reinforcement  agents  [15-17],  Though 
the  mechanical  compressive  performance  of  CNTs  is 
expected  to  contribute  significantly  to  the  performance  of 
these  applications,  relatively  little  is  known  about  the 
collective  mechanical  behavior  of  CNT  arrays  and  how  their 
morphology  may  influence  mechanical  response.  Although 
individual  CNTs  have  demonstrated  experimentally  mea¬ 
sured  tensile  elastic  modulus  values  on  the  order  of  1  TPa 
[18,19],  experimentally  obtained  compressive  modulus 
values  for  vertically  aligned  CNT  arrays  are  often  on  the 
order  of  10  MPa  -  a  difference  of  five  orders  of  magnitude. 


Because  of  the  inherent  waviness  and  strong  van  der  Waals 
interactions  between  neighboring  CNTs  [20],  classical  Euler- 
Bernoulli  beam  theory  alone  is  insufficient  for  capturing 
the  complex  behavior  of  the  networked  interactions  present 
within  CNT  arrays  [21].  Current  mechanical  models  of  CNT 
lack  the  incorporation  of  the  complexities  of  CNT  morphol¬ 
ogy,  and  additional  insight  into  the  compressive  behavior 
of  these  networks  is  required. 

Nanoindentation  is  a  popular  means  for  material  property 
evaluation  for  thin  films  and  coatings.  The  method  employs 
high  precision  instrumentation  to  monitor  force  and  dis¬ 
placement  of  an  indenter  tip  during  compressive  loading 
and  unloading  of  a  sample.  Elastic  properties  are  typically 
obtained  from  the  force-displacement  slope  (stiffness)  during 
initial  withdrawal  of  the  tip,  where  forces  arise  due  to 
elastic  recovery  alone.  Discrete  evaluations  of  this  type  are 
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considered  valid  only  at  the  depth  of  initial  tip  withdrawal.  To 
evaluate  elastic  properties  as  a  continuous  function  of  pene¬ 
tration  depth,  the  continuous  stiffness  measurement  (CSM) 
technique  introduces  relatively  high  frequency  loading  and 
unloading  cycles  by  imparting  a  small  sinusoidal  displace¬ 
ment  at  the  indenter  tip  superimposed  upon  the  larger  steady 
tip  displacement  rate.  In  such  a  way,  stiffness  and  elastic 
properties  may  be  evaluated  at  each  data  collection  interval. 
Multiple  thorough  nanoindentation  reviews  exist  in  the  liter¬ 
ature  detailing  the  theoretical  basis  of  the  technique,  effects 
of  tip  geometry,  and  applications  [22-25]. 

Several  researchers  have  employed  nanoindentation  and 
similar  instrumented  indentation  methods  as  means  to  eval¬ 
uate  the  compressive  mechanical  properties  of  CNT  arrays 
[26-34],  the  results  of  which  are  summarized  in  Table  1. 
Tested  CNT  morphologies  include  short,  non-interacting  ar¬ 
rays  [26,27],  typical  CVD-grown  arrays  with  heights  from 
15  |jm  to  over  1  mm  [28-33],  and  densely  populated  and  highly 
aligned  CNT  brushes  [34].  Data  collection  techniques  and 
hardware  inconsistencies  make  bridging  data  sets  difficult, 
and  as  a  result  no  clear  trends  may  be  established  relating 
the  mechanical  properties  to  CNT  array  parameters  such  as 
length,  diameter,  and  morphology  based  upon  the  available 
data.  Compressive  elastic  properties  within  the  literature 
are  evaluated  most  commonly  using  three  techniques.  These 
include  evaluations  based  upon  the  initial  stress-strain  load¬ 
ing  slope,  tangential  unloading  stiffness,  or  the  CSM  tech¬ 
nique.  Because  of  the  elastic-plastic  behavior  of  most  CNT 
arrays,  the  loading  slope  may  be  significantly  lower  than  that 
of  the  unloading  slope,  directly  leading  to  differences  in  cal¬ 
culated  modulus  results.  To  further  confound  data,  indenta¬ 
tion  tip  geometries  vary  between  experiments  and  include 
spherical,  pyramidal,  Berkovich,  rectangular  struts,  and  flat 
punch  geometries.  The  influence  of  tip  geometry  relative  to 
the  elastic  behavior  of  CNT  arrays  is  not  well  studied,  though 
unintended  transverse  forces  will  undoubtedly  be  imparted  to 
the  CNT  tips  utilizing  angular  tips.  Nevertheless,  most  evalu¬ 
ations  report  an  array  elastic  compressive  modulus  of  less 
than  100  MPa.  A  thorough  parametric  study  of  CNT  array  mor¬ 
phology  may  help  to  decouple  trends  associated  with  experi¬ 
mental  technique  from  those  generated  by  the  CNT  arrays 
themselves. 

The  mechanical  yield  of  indented  CNT  arrays,  defined  as  a 
reduction  in  the  stress-strain  slope  after  initial  linear  loading, 
has  likewise  been  observed  during  indentation  [26-29,34],  For 
those  experiments  producing  yielding  behavior,  typical  buck¬ 
ling  yield  strains  range  from  0.01  to  0.6,  as  observed  in  Table  1. 
Others  have  evaluated  only  the  top  several  microns  of  an  ar¬ 
ray,  a  penetration  depth  that  presumably  precedes  the  onset 
of  yielding  [30-33].  Unfortunately,  scanning  electronic  micro¬ 
scope  (SEM)  evaluation  of  yielded  arrays  after  indentation  is 
often  not  reported,  and  the  morphology  change  of  CNTs 
resulting  from  the  yielding  behavior  is  unclear.  Zbib  et  al.  re¬ 
ported  micron-scale  wrinkles  localized  at  the  base  of  defined 
CNT  pillars  arose  uniaxial  compression  [28],  Significant  lat¬ 
eral  translation  at  the  base  of  the  pillars  occurred  during  com¬ 
pression,  and  it  is  unclear  how  the  translation  influenced  the 
wrinkle  formation.  As  with  the  elastic  mechanical  property 
evaluation,  consistent  trends  with  respect  to  mechanical 
yielding  from  existing  literature  are  not  apparent,  and  the 


diversity  of  indentation  techniques  further  confounds  the 
ability  to  compare  results.  As  a  consequence,  the  physical 
mechanisms  governing  the  yielding  behavior  of  indented 
CNT  arrays  remain  poorly  understood,  and  the  specific  role 
of  CNT  array  height  over  several  orders  of  magnitude  with 
respect  to  mechanical  yielding  behavior  has  not  been 
explored. 

Freestanding  CNT  blocks,  removed  from  a  native  growth 
substrate,  have  similarly  been  examined  in  uniaxial  compres¬ 
sion.  Reported  elastic  compressive  modulus  values  vary  from 
0.9  to  20.8  MPa  -  a  much  narrower  distribution  than  those 
reported  by  localized  indentation  of  substrate-adhered  CNT 
arrays.  Typically,  compression  of  a  freestanding  CNT  block 
takes  place  between  parallel  compression  cylinders,  resulting 
in  uniaxial  compression  over  the  entire  cross  section  of  the 
block,  in  contrast  to  the  localized  compression  imparted  by 
nanoindentation.  Freestanding  CNT  blocks  may  be  generated 
from  direct  liberation  from  a  native  growth  substrate  [35-37] 
or  from  compression  molding  of  loose  CNTs  from  solution 
into  freestanding  columns  [38].  Directly-liberated  films  dem¬ 
onstrate  remarkable  elastic  recovery,  even  after  rigorous 
mechanical  cycling.  An  860  pm  liberated  array  was  shown  to 
recover  its  full  height  after  85%  compressive  strain,  with 
nearly  full  recovery  observed  for  hundreds  of  mechanical  cy¬ 
cles  [35],  After  10,000  compressive  cycles  of  85%  strain  ampli¬ 
tude,  the  array  was  lost  only  20%  of  its  initial  height.  Gradual 
height  loss  within  the  CNT  blocks  was  attributed  to  the  for¬ 
mation  and  evolution  of  highly  coordinated  parallel  buckling 
folds  observed  at  the  bottom  of  the  array  and  originating  after 
200  cycles.  In  a  separate  experiment,  500,000  cycles  of  15% 
strain  amplitude  resulted  in  a  height  reduction  of  less  than 
5%  for  a  2.2  mm  liberated  array  [36],  Coordinated  buckles  were 
also  observed  after  compression  of  a  3.4  mm  liberated  CNT 
block  [37].  Although  the  CNT  buckles  were  on  the  order  of 
10  pm  in  wavelength,  their  formation  was  not  readily 
detected  within  the  resulting  stress-strain  behavior.  In 
contrast,  centimeter-scale  pillars  cast  from  CNT  solutions 
recovered  to  only  40%  strain  after  a  single  loading  to  85% 
strain,  with  decreased  recovery  observed  at  increased  cycles 
[38],  These  drastic  differences  in  compressive  response  illus¬ 
trate  a  strong  coupling  between  CNT  array  morphology  that 
is  currently  unexplained. 

Within  this  paper,  we  report  the  uniaxial  nanoindentation 
of  CNT  arrays  of  similar  morphology  with  heights  differing  by 
more  than  an  order  of  magnitude.  By  utilizing  a  flat  punch 
nanoindentation  tip,  hundreds  of  microns  of  penetration 
depth  were  achievable  without  imparting  unintended  lateral 
forces.  Nanoindentation  of  substrate-adhered  CNT  arrays  to 
these  depths  has  not  been  previously  reported.  The  use  of 
the  CSM  technique  allowed  for  continuous  elastic  property 
evaluation  throughout  the  indentation  depth.  Lateral  force 
measurements  were  also  recorded  as  a  function  of  penetra¬ 
tion  depth,  yielding  a  three-dimensional  profile  of  forces  aris¬ 
ing  during  indentation.  Metrics  of  interest  for  the  evaluation 
include  compressive  forces,  lateral  forces,  elastic  stiffness, 
and  changes  in  CNT  morphology  observed  using  SEM 
analysis.  The  resulting  mechanical  response  of  the  arrays  is 
compared  to  that  of  open-cell  foam,  and  the  resulting  defor¬ 
mation  modes  are  described  utilizing  concepts  of  cellular 
materials. 


Table  1  -  Reported  compressive  CNT  mechanical  properties. 

Support  substrate 

Indenter  tip/compressive 

CNT  height 

CNT  diameter 

Onset  of  buckling 

Elastic  compressive 

Modulus  computation 

technique 

strain 

modulus 

method 

Anodic  Alumina  [24] 

2  pm  Flat  punch 

50,  100  nm 

25  nm 

0.2,  0.6 

N/A 

N/A 

Anodic  Alumina  [25] 

2  pm  Flat  punch/100  pm  diameter 

50,  100  nm 

50  nm 

0.2,  0.6  flat  punch/N/A, 

N/A 

N/A 

sphere/100  nm  Berkovich 

0.4  sphere/0.05-0.08, 

0.15  Berkovich 

Silicon  [26] 

Berkovich  -  1000  nm  diameter  tip 

60  pm 

N/A 

0.03 

14.9  ±  5.7  MPa 

Unloading  stiffness 

Adhered  to  solid  [27] 

1  mm  Diameter  probe 

>1  mm 

43  nm 

0.1 

27.6  MPaa 

Loading  stiffness 

Silicon  [28] 

Quartz  strut  -  2  x  0.5  mm 
cross  section 

15,  40,  500  pm 

20-30  nm 

N/A 

0.22-0.25  MPa 

Loading  stiffness 

Polymer  [29] 

30  pm  Flat  punch 

>600  pm 

15  nm 

N/A 

N/A 

N/A 

Crystalline  [30] 

3-Sided  pyramid 

570-1100  nm 

55-104  nm 

N/A 

0.9-1.23  Tpa 

Bending  stiffness 

Silicon  [31] 

Berkovich 

N/A 

N/A 

N/A 

40-600  MPa 

Continuous  stiffness 

SiC  [32] 

2,  10,  27  pm  Diameter  spheres 

0.2  and  1.3  pm 

1-3  nm 

0.08,  0.03,  0.01 

18  GPa 

Loading  stiffness 

Freestanding  film  [33] 

Parallel  plates 

3.4  mm 

53  nm 

N/A 

0.9  MPac/2.2  MPad 

Loading  stiffness 

Freestanding  film  [34] 

Parallel  plates 

2.2  mm 

50  nm 

0.01b 

2.1  MPac/20.8  MPad 

Loading  stiffness 

Freestanding  film  [35] 

Parallel  plates 

860  pm 

40  nm 

0.22 

12  MPa 

Loading  stiffness 

Freestanding  pillar  [36] 

Parallel  plates 

>1  cm 

N/A 

0.15 

10-12  MPa 

Loading  stiffness 

a  Value  removes  3%  area  normalization. 

b  First  mechanical  cycle. 

c  Measured  at  initial  loading  condition. 
d  Measured  after  plateau  region. 
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Raman  Shift  (cm  ) 


Fig.  1  -  (a)  SEM  micrograph  of  100  |im  diameter  flat  punch  indenter.  (b)  Raman  spectra  obtained  from  indented  CNT  arrays.  The 
inset  numbers  indicate  the  CNT  array  heights. 


2.  Experimental  methods 

A  nominally  100  pm  diameter  diamond  flat  punch  (Micro  Star 
Technologies)  was  selected  as  the  nanoindentation  tip  and  is 
shown  in  Fig.  la.  The  relatively  large  surface  area  (9001  pm2) 
allows  for  sensitive  detection  of  the  CNT  top  surface,  a  nearly 
uniform  contact  area  with  penetration,  and  uniaxial  com¬ 
pression.  The  flat  punch  features  a  taper  of  less  than  5°  to  a 
depth  of  approximately  300  pm  from  the  indenter  face.  This 
geometry  allows  for  indentation  depths  spanning  hundreds 
of  microns  with  little  geometrical  influence. 

Multiwalled  CNT  arrays  were  synthesized  on  15  x  15  mm 
Si02/Si  substrates  using  two  procedures.  A  1200  pm  array 
was  synthesized  utilizing  a  water-assisted  CVD  technique 
incorporating  ethylene  as  a  carbon  feedstock  with  hydrogen 
and  helium  carrier  gasses  [39].  This  addition  of  water  vapor 
has  been  shown  to  reduce  amorphous  carbon  on  catalyst  par¬ 
ticles,  extend  catalyst  lifetime,  and  produce  CNT  arrays  in  ex¬ 
cess  of  1  mm  [40],  For  this  procedure  a  1  nm  thick  Fe  catalyst 
layer  was  deposited  atop  a  10  nm  thick  Al203  buffer  film. 
Water  vapor  was  introduced  at  concentrations  ranging  from 
200  to  400  ppm.  Synthesis  proceeded  at  1  atm  and  760  °C  for 
ten  hours  to  generate  an  array  height  of  1200  pm.  CNTs  pro¬ 
duced  by  this  procedure  exhibited  a  diameter  distribution 
between  10-15  nm  and  3-8  walls,  as  observed  using  transmis¬ 
sion  electron  microscope  (TEM)  analysis  (see  Supplementary 
data).  The  remaining  arrays  were  produced  by  utilizing  cata¬ 
lyst  film  with  3  nm  Fe  atop  10  nm  of  Al  in  the  absence  of  water 
[41].  This  procedure  yields  CNT  arrays  of  various  heights 
depending  upon  the  synthesis  conditions.  CNT  synthesis 
was  achieved  at  a  growth  temperature  of  750  °C  by  flowing  ar¬ 
gon,  hydrogen,  and  acetylene  at  10-100  Torn  CNT  arrays  were 
synthesized  to  heights  of  35,  190,  300,  and  650  pm  with  this 
procedure.  TEM  analysis  revealed  that  CNTs  consisted  of  2-3 
walls  (see  Supplementary  data),  with  outer  diameters  be¬ 
tween  10  and  20  nm.  Raman  spectra  of  the  CNT  arrays,  dis¬ 
played  in  Fig.  lb,  were  obtained  by  a  Renishaw  spectrometer 
operating  at  a  laser  wavelength  of  532  nm.  The  ratio  of  the 
D-band  (1340  cm-1)  to  G-band  (1580  cmT1)  for  all  samples 
gives  an  indication  that  the  arrays  were  of  comparable  crys¬ 
talline  quality  and  composition  as  measured  in  the  top  and 


off-axis  direction  with  an  unpolarized  source  (see  Supporting 
information).  Typical  ratios  ranged  from  0.68  to  0.83,  with  the 
1.2  mm  array  exhibiting  the  greatest  D:G  ratio. 

Instrumented  indentation  was  performed  using  an  MTS 
Nanoindenter  XP  system  equipped  with  lateral  force  and 
CSM  capabilities.  The  CSM  parameters  utilized  in  this  study 
include  a  2  nm  amplitude  and  50  Hz  frequency,  consistent 
with  recommendations  detailed  in  the  MTS  machine  manual 
and  similar  to  conditions  used  by  others  [34].  Surface  contact 
was  defined  as  the  position  at  which  the  compressive  force- 
displacement  slope  exceeded  25  N/m  during  tip  approach. 
Stress  was  controlled  such  that  the  loading  rate  divided  by 
the  applied  load  was  held  at  a  constant  of  0.01  s-1.  Such  a 
loading  scheme  is  approximately  equivalent  to  that  of  a  con¬ 
stant  indentation  strain  rate  with  materials  exhibiting  con¬ 
stant  hardness  with  depth  [42],  A  typical  duration  for  a 
single  indent  was  approximately  3-5  min.  Data  was  collected 
by  indenting  no  fewer  than  20  locations  within  a  given  array 
with  grid-like  spacing  of  200-500  pm  between  neighboring  in¬ 
dents.  Spacing  between  indents  had  no  impact  on  mechani¬ 
cal  response  within  this  separation  range.  Indents  were 
ceased  upon  either  reaching  a  maximum  displacement  of 
200  pm  for  arrays  having  a  height  of  greater  than  300  pm,  or 
upon  reaching  a  defined  compressive  load  of  200  mN  for 
shorter  arrays.  The  indenter  tip  was  cleaned  between  tests 
by  subjecting  the  tip  to  a  series  of  lateral  scratch  tests  on  sil¬ 
ica,  which  removed  any  CNT  debris  that  may  adhere  to  the  tip 
during  compression. 


3.  Results  and  discussion 

To  validate  the  efficacy  of  the  CSM  technique  with  CNT  arrays, 
the  elastic  stiffness  obtained  using  CSM  was  compared  to  that 
obtained  by  measuring  the  tangential  unloading  stiffness  at 
discrete  depths.  A  300  pm  CNT  array  was  chosen  for  this  eval¬ 
uation,  with  10  locations  sampled  for  each  method  to  a  depth 
of  200  pm.  The  discrete  tangential  stiffness  measurements 
from  the  force  unloading  curves  were  obtained  at  approxi¬ 
mately  every  20  pm  of  depth.  Stiffness  analyses  were  com¬ 
puted  by  fitting  a  linear  slope  from  the  beginning  of  the 
unloading  segment  to  the  point  of  40%  unloading  for  both 
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Fig.  2  -  (a)  Stiffness  measured  from  300  pm  CNT  array  using  CSM  (lines)  and  tangential  unloading  (filled  circles)  techniques,  (b) 
Load  vs.  displacement  profile  for  a  typical  discrete  tangential  stiffness  measurement,  (c)  Load-displacement  behavior  of  a 
300  |im  CNT  array  for  100  mechanical  cycles  and  (d)  the  resulting  elastic  stiffness  as  a  function  of  mechanical  cycle. 


methods.  Because  the  unloading  slope  is  not  strictly  linear, 
altering  the  unloading  point  at  which  stiffness  measurements 
were  calculated  resulted  in  minor  changes  in  measured  stiff¬ 
ness  values,  with  the  changes  remaining  relatively  consistent 
between  both  methods.  Prior  to  alleviation  of  applied  load, 
the  load  was  held  constant  for  5  s  to  minimize  potential  creep 
effects.  The  resulting  stiffness  profiles  are  shown  in  Fig.  2a 
and  show  strong  agreement  with  data  obtained  using  the  dis¬ 
crete  tangential  stiffness  method.  CSM  data  is  displayed  as 
continuous  line  profiles,  while  discrete  tangential  stiffness 
data  are  displayed  closed  circles.  Fig.  2b  shows  significant 
hysteresis  between  loading  and  unloading  cycles,  consistent 
with  observations  of  others  [43].  Though  the  CNT  arrays  re¬ 
cover  greater  than  1  pm  of  compression  during  discrete  tan¬ 
gential  unloading  cycles,  compared  to  less  than  2  nm 
recovery  during  CSM  during  stiffness  evaluation,  a  strong 
agreement  exists  between  the  two  measurement  techniques. 
As  a  result,  CSM  measurements  served  as  the  basis  for  com¬ 
pressive  mechanical  property  evaluation  for  the  remainder 
of  experimental  analyses,  unless  otherwise  stated. 

To  examine  the  repeatability  of  elastic  stiffness  with 
mechanical  cycling,  a  300  pm  CNT  array  was  loaded  to  a  pen¬ 
etration  depth  of  200  pm,  at  which  point  the  load  was  reduced 
to  5%  of  the  maximum  load.  The  load  was  then  reapplied  un¬ 
til  the  200  pm  depth  was  again  achieved,  and  the  process  was 


repeated  for  100  total  cycles.  A  nearly  constant  elastic  stiff¬ 
ness  is  maintained  with  respect  to  mechanical  cycling.  Com¬ 
pressive  loads  were  maintained  at  the  maximum  and 
minimum  values  for  5  s  at  the  end  of  the  unloading  and  load¬ 
ing  cycles  to  accommodate  creep  processes.  Fig.  2c  and  d 
shows  the  representative  mechanical  cycling  behavior  and 
the  resulting  tangential  unloading  stiffness  calculated  at 
40%  unload.  The  hysteresis  loop  between  loading  and  unload¬ 
ing  cycles  evolves  significantly  between  first  and  second  cy¬ 
cles,  and  gradually  reaches  a  repeatable  state  after 
approximately  50  cycles.  Similar  hysteretic  behavior  has  been 
observed  by  others  [34,38]  and  has  been  attributed  to  energy 
dissipation  from  the  formation  and  breaking  of  neighbor-to- 
neighbor  bonds  or  viscous  loss  from  air  flow  induced  during 
compression.  Note  that  the  array  recovered  only  8  pm  (4%) 
of  compression  after  the  initial  compression  cycle,  with 
increasingly  less  recovery  observed  at  advanced  cycles.  This 
corresponds  to  than  96%  retained  plastic  deformation.  As 
demonstrated  in  Fig.  2d,  the  stiffness  increases  by  less  than 
15%  between  the  first  and  hundredth  cycle,  with  a  repeatable 
stiffness  maintained  after  approximately  50  cycles.  The 
unloading  stiffness  evolution  is  related  to  the  strain  recovery 
behavior  in  that  progressively  less  strain  recovery  necessi¬ 
tates  an  increased  stiffness  at  the  point  of  unloading.  After 
a  repeatable  strain  recovery  is  established,  a  constant 
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stiffness  value  of  approximately  1680  N/m  is  likewise  ob¬ 
served.  This  behavior  demonstrates  the  potential  to  tune  an 
array’s  density  and  slightly  alter  the  mechanical  stiffness. 

3.1.  Stress-strain  behavior 

Compressive  stress  is  defined  as  the  compressive  force  di¬ 
vided  by  the  indenter  tip  face  area  of  9001  pm2.  The  flat  punch 
geometry  used  in  this  study  simplifies  the  computation  of 
mechanical  stress  compared  to  indenters  having  spherical 
or  angular  geometries,  as  the  cross-sectional  area  of  these 
types  of  indenters  is  a  strong  function  of  indentation  depth. 
Similarly,  strain  was  computed  by  dividing  the  compressive 
displacement  by  the  initial  CNT  array  height.  These  straight¬ 
forward  and  intuitive  definitions  of  stress  and  strain  are  an 
additional  advantage  of  utilizing  a  flat  punch  tip  geometry. 
Although  proper  alignment  between  the  face  of  a  flat  tip  in¬ 
denter  and  the  sample  is  a  concern,  alignment  co-planarity 
differences  were  measured  at  less  than  1  angular  degree.  Be¬ 
cause  of  the  significant  depth  of  penetration  for  this  study, 
and  in  light  of  potential  inherent  height  variation  within  the 
free  tip  area,  any  misalignment  between  the  indenter  tip 
and  free  surface  of  the  arrays  is  considered  to  be  small  and 
isolated  to  only  the  top  2  pm  of  contact. 

The  collective  compressive  response  of  the  arrays,  ob¬ 
served  in  Fig.  3,  consists  of  three  distinct  stress-strain  regions 
and  is  consistent  with  that  of  an  open-cell  foam,  as  was  first 
reported  by  Cao  et  al.  [35].  The  initial  region  features  a  linear 
or  near-linear  region  at  low  stress,  followed  by  an  elongated 
cellular  collapse  plateau,  and  a  final  densification  region  in 
which  stress  rapidly  increases.  The  plateau  region  for  tradi¬ 
tional  foam  may  consist  of  recoverable  elastic  buckling,  plas¬ 
tic  yielding,  or  brittle  crushing,  depending  on  the  foam 
composition  and  mechanical  properties  of  the  cellular  strut 
material.  All  three  mechanisms  result  in  a  similar  stress- 
strain  behavior.  Mechanical  behavior  in  each  region  is  dic¬ 
tated  in  large  part  by  the  relative  density  of  the  strut  constit¬ 
uent  material.  The  relative  density  of  cellular  strut  material  is 
strongly  related  to  the  thickness  and  geometry  of  the  struts  by 
the  following  relationship 


Strain  (pm/pm) 


Fig.  3  -  Average  stress-strain  behavior  of  the  indented  CNT 
arrays.  The  inset  numbers  indicate  the  CNT  array  heights. 


where  p  is  the  mass  density,  t  is  the  open  cell  strut  thickness,  1 
is  the  characteristic  length  of  a  strut  in  the  direction  normal 
to  length,  the  subscript  s  refers  to  the  properties  of  the  strut 
solid  material,  and  the  superscript  p*  refers  to  the  observed 
properties  of  the  foam  material.  Proportionality  may  be  deter¬ 
mined  from  the  cellular  geometry. 

The  onset  of  elastic  and  plastic  yield  of  cellular  solids  are 
strongly  dependent  on  relative  density  and  follow  the 
relations 


for  elastic  foam  material  and 


for  plastically  yielding  foam,  where  <r*,  is  the  elastic  collapse 
stress,  E  is  the  elastic  modulus,  <7*  represents  the  plastic  col¬ 
lapse  stress,  <7ys  is  the  yield  strength  of  the  cellular  strut  mate¬ 
rial,  and  Ci  and  C2  are  constants.  For  low  density  elastic  foam 

<  0/3),  a  constant  of  proportionality  (Ca)  of  0.05  empirically 
fits  most  materials.  Similarly,  the  predicted  collapse  strain  for 
elastic  foam  is  simply  equal  to  the  constant  of  proportionality 
of  0.05.  A  value  of  0.03  has  been  experimentally  determined 
for  C2,  resulting  in  a  predicted  plastic  collapse  strain  similar 
to  that  of  elastic  foam  at  moderate  densities  and  slightly 
higher  for  lower  density  plastic  foam.  Similarly,  the  onset  of 
densification  strain  in  an  open-cell  foam  decreases  as  a  func¬ 
tion  of  relative  density. 

The  average  compressive  behavior  of  CNT  arrays  having 
various  heights,  displayed  in  Fig.  3,  is  consistent  with  that 
of  an  open-cell  foam  material  of  various  relative  densities. 
The  35  pm  array  behaves  similarly  to  low  density  foam,  while 
the  1200  pm  array  exhibits  the  behavior  of  higher  density 
foam.  Intermediate  array  heights  also  follow  the  general  rela¬ 
tionship  that  greater  CNT  heights  exhibit  behavior  of  in¬ 
creased  foam  density.  The  average  yield  stress  ranges  from 
approximately  0.01  MPa  for  the  35  pm  array  to  2.9  MPa  for 
the  1200  pm  array.  Yield  strain  occurs  at  approximately  0.03 
for  the  1200  pm  array  to  0.12  for  the  35  pm  array,  with  a  trend 
similar  to  that  anticipated  for  plastically  yielding  foam  [44], 
The  onset  of  densification  strain  decreases  with  array  height, 
consistent  with  open-cell  foam  behavior  at  various  densities. 

Closer  examination  of  individual  stress-strain  curves, 
shown  in  Fig.  4,  yields  additional  insight  into  the  compressive 
behavior  of  each  array  and  demonstrates  the  run-to-run  var¬ 
iability  observed  for  each  sample.  The  35  pm  array  exhibited 
the  greatest  run-to-run  variability  with  respect  to  the  yield 
strain  and  densification  strain,  perhaps  due  in  part  to  diffi¬ 
culty  associated  with  accurate  force-based  detection  of  the  ar¬ 
ray  surface.  It  is  noteworthy  to  emphasize  that  the  35  pm 
array  demonstrated  significantly  lower  initial  loading  slope 
relative  to  the  other  arrays.  Interestingly,  all  samples  other 
than  the  35  pm  array  exhibit  a  periodic  stress-strain  oscilla¬ 
tion  signature  within  the  plateau  region  with  characteristic 
wavelengths  between  approximately  3-6  pm.  To  investigate 
the  physical  mechanism  giving  rise  to  the  stress-strain  oscil¬ 
lations,  each  CNT  array  was  indented  along  a  free  edge  of  the 
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Fig.  4  -  Stress-strain  relationship  for  indented  (a)  35  jim,  (b)  190  (im,  (c)  300  (im,  (d)  650  (im,  and  (e)  1200  fim  arrays. 


Fig.  5  -  Characteristic  buckling  behavior  of  (a)  35  jim,  (b)  190  nm,  (c)  650  jim,  and  (d)  1200  (im  CNT  arrays.  Scale  bars  of  inset 
images  =  10  |im. 

substrate  support  to  a  depth  exceeding  the  onset  of  the  yield  samples  is  consistent  with  that  obtained  within  the  bulk  of 
strain.  The  stress-strain  behavior  at  the  free  edges  of  the  the  sample.  This  elicits  confidence  that  the  observed  yielding 
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Fig.  6  -  Buckling  characteristics  of  650  |im  CNT  array,  (a)  SEM  micrograph  of  buckled  region  of  75  nm  indent  with  inset 
showing  undisturbed  CNTs  directly  under  buckling  front,  (b)  Corresponding  instantaneous  loading  slope  for  75  |im  indent. 


behavior  is  an  intrinsic  characteristic  of  the  CNTs,  indepen¬ 
dent  from  neighboring  interactions.  Fig.  5  shows  cross-sec¬ 
tional  SEM  images  of  the  indented  edges  which  reveal  clear 
evidence  of  plastic  hinge  and  buckle  formation  for  each  array, 
consistent  with  plastic  yielding  of  cellular  struts  within  a 
foam  material. 

The  35  pm  array  displayed  an  out  of  plane,  S-shaped  bend¬ 
ing  behavior  over  the  height  of  the  array.  The  buckle  is  pre¬ 
sumed  to  originate  during  the  plateau  region  of  the  stress- 
strain  curve.  Arrays  having  heights  of  190,  300,  and  650  fim 
display  a  different  behavior,  with  similar  characteristics. 
These  arrays  exhibited  localized,  accordion-like  collective 
buckling  and  folding  originating  at  the  top  surface  of  the  in¬ 
dented  array.  Increasing  indentation  depths  resulted  in  the 
formation  of  a  greater  number  of  buckling  folds.  No  evidence 
of  compression  existed  within  the  array  region  directly  be¬ 
neath  the  buckling  front  for  these  arrays,  indicating  highly 
localized  compression  within  the  buckling  region.  However, 
evidence  of  a  single  buckle  at  the  bottom  of  the  array,  near 
the  growth  substrate,  was  observed  after  indentation  of  the 
650  |im  array.  The  bottom  buckle  is  unique  in  appearance 
compared  to  the  buckles  appearing  near  the  top  surface  and 
appears  to  have  recovered  a  significant  portion  of  deflection. 
The  origination  of  coordinated  buckles  near  the  arrays’  top 
surfaces  contradicts  results  obtained  from  the  compression 
of  liberated  millimeter-length  CNT  arrays,  in  which  coordi¬ 
nated  buckling  folds  were  observed  strictly  on  the  bottom  side 
of  the  CNT  film  [35,37].  Additionally,  the  deformation  ob¬ 
served  herein  by  nanoindentation  is  largely  plastic  in  nature, 
with  less  than  5%  strain  recovery  achieved  within  the  plateau 
region  after  one  cycle,  while  uniaxial  compression  of  liber¬ 
ated  arrays  is  largely  elastic  in  nature  [35].  Furthermore, 
SEM  examination  of  the  1200  pm  array  reveals  clearly  plastic 
deformation,  with  much  of  the  indent  strain  unrecovered.  De¬ 
spite  the  occurrence  of  periodic  stress-strain  oscillations  dur¬ 
ing  indentation,  the  buckling  behavior  lacks  the  organization 
and  coordination  exhibited  by  other  arrays.  However,  isolated 
small-amplitude  ripples  spanning  a  portion  of  the  projected 
indent  area  are  observed  at  various  depths  beneath  the  sur¬ 
face.  A  lack  of  large-scale  buckling  with  largely  plastic  defor¬ 
mation  may  indicate  densification  of  CNTs  throughout  a 
greater  depth  of  the  array.  Indeed,  Fig.  5d  shows  that  small- 


scale  buckles  are  not  localized  at  the  top  of  the  indented  sur¬ 
face,  as  is  the  case  with  other  arrays.  An  additional  large 
wavelength  buckle  is  observed  at  the  bottom  of  the  array  (di¬ 
rectly  above  the  substrate  surface).  The  observed  buckle  is 
similar  to  that  observed  after  indentation  of  the  650  pm  array, 
but  larger  in  wavelength.  In  the  absence  of  in  situ  observation, 
it  is  difficult  to  determine  how  buckling,  densification,  and 
deformation  evolved  within  this  array.  However,  the  compres¬ 
sive  behavior  of  this  array  is  unique  amongst  the  others 
tested  and  may  originate  in  part  from  the  differentiated  syn¬ 
thesis  method  and  differences  in  wall  count.  Lower  magnifi¬ 
cation  images  of  the  650  and  1200  pm  indented  arrays  may 
be  found  in  Supplementary  data  section. 

To  demonstrate  the  relationship  between  the  force-dis¬ 
placement  oscillations  and  the  formation  of  mechanical 
buckles,  the  number  of  stress-strain  oscillations  was  com¬ 
pared  to  the  number  of  observed  mechanical  folds.  For  this 
examination,  the  650  pm  array  was  indented  to  a  depth  of 
75  pm  (0.12  strain)  along  a  free  edge  of  the  sample  using  stan¬ 
dard  parameters.  SEM  analysis,  shown  in  Fig.  6,  displays  se¬ 
ven  visible  folds  after  indentation  (the  top  region  above  the 
first  trough  shows  no  signs  of  having  buckled).  To  more 
clearly  analyze  the  number  of  force  oscillations  recorded  dur¬ 
ing  the  nanoindentation,  the  instantaneous  loading  slope, 
displayed  in  Fig.  6b,  was  evaluated  to  magnify  the  amplitude 
of  oscillations  over  that  observed  by  force-displacement.  The 
plot  reveals  14  force  oscillations,  corresponding  to  double  the 
number  of  folds  observed  by  SEM  analysis.  Based  on  this 
information,  we  surmise  that  a  complimentary  set  of  buckles 
exists  on  the  side  opposite  the  free  edge  observed  in  Fig.  6a, 
with  each  complimentary  buckle  exhibiting  a  fold  out  of 
phase  with  those  visible  at  the  free  edge  to  form  an  “S- 
shaped”  fold. 

The  formation  of  buckles  within  the  stress  collapse  pla¬ 
teau  is  conceptually  akin  to  the  formation  of  plastic  hinges 
exhibited  within  the  cells  of  traditional  open-cell  foam  under¬ 
going  plastic  collapse  [44],  Hinges  form  when  the  localized 
stress  on  cell  struts  is  greater  than  the  strut  material’s  com¬ 
pressive  yield  strength,  leading  to  plastic  cellular  collapse. 
Plastic  buckling  would,  then,  presumably  occur  first  within 
those  cells  most  prone  to  yielding.  Combining  Eqs.  (1)  and 
(3)  reveals  that  low  plastic  collapse  stress  results  from  low 
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Fig.  7  -  Schematic  representation  of  CNT  alignment  as  a  function  of  depth  and  the  resulting  compressive  buckling  formation 
as  a  function  of  array  height. 


cellular  density,  which  is  itself  highly  sensitive  to  cellular 
geometry,  as  illustrated  by  the  relationship 


Assuming  consistent  cellular  shape,  Eq.  (4)  suggest  that  cells 
having  a  broad  span  perpendicular  to  axial  growth  direction 
(corresponding  to  low  density  and  perpendicular  alignment) 
would  exhibit  the  lowest  yield  stress.  To  that  end,  Bedewy 
et  al.  [45]  recently  spatially  resolved  the  vertical  alignment 
of  CNTs  within  arrays  utilizing  synchrotron  X-ray  scattering. 
Based  on  their  analysis  of  the  Hermans  orientation  parame¬ 
ter,  vertical  CNT  alignment  and  density  were  observed  to  in¬ 
crease  with  depth  from  the  top  surface  of  arrays  until  an 
inflection  point  at  a  depth  of  a  few  hundred  micrometers 
was  reached,  at  which  point  alignment  decreased  monotoni- 
cally  with  depth.  The  density  and  vertical  alignment  at 
approximately  500  pm  depth  had  reduced  sufficiently  to  be 
equivalent  to  that  observed  on  the  surface,  and  decreased  fur¬ 
ther  with  depth  beyond  500  pm,  reaching  a  minima  near  the 
growth  substrate.  CNT  arrays  whose  growth  was  terminated 
prior  to  attaining  500  pm  height  exhibited  a  similar  Hermans 
orientation  parameter  profile,  featuring  a  maxima  near 
200  pm  of  depth  and  minima  near  the  top  surface.  Translating 
these  results  to  the  CNT  arrays  studied  herein,  we  may  de¬ 
duce  that  cellular  elongation  and  density  (and  related  plastic 
collapse  stress)  is  lowest  at  the  top  surface  for  all  arrays  other 
than  the  650  and  1200  pm  arrays,  where  the  minima  resides 
near  the  bottom  surface.  Raman  spectroscopy  analysis  has 
likewise  revealed  that  structural  defects  in  super-long  CNT  ar¬ 
rays,  similar  in  structure  to  the  1200  pm  array,  increase  as  a 
function  of  distance  from  the  array  top  surface.  Arguments 
based  upon  cellular  alignment  are  also  consistent  with  the 
buckling  behavior  observed  for  liberated  arrays,  which  buckle 
solely  on  the  bottom  surface.  As  illustrated  in  Fig.  7,  localiza¬ 


tion  of  coordinated  buckling  near  the  top  surface  for  arrays  of 
moderate  heights  and  additional  buckling  near  the  substrate 
of  the  650  and  1200  pm  arrays  may  be  explained  due  to  these 
factors. 

3.2.  Lateral  force  measurement 

Lateral  force  measurement  may  provide  additional  insight 
both  into  the  magnitude  and  relative  direction  of  cellular  col¬ 
lapse  during  the  buckling  events.  Fig.  8  shows  representative 
lateral  and  normal  force  data  obtained  from  the  650  to 
1200  pm  arrays.  Lateral  force  oscillations  arise  at  a  strain  con¬ 
sistent  with  the  onset  of  normal  force  oscillations,  indicating 
lateral  motion  is  generated  during  buckle  formation. 
Although  the  absolute  value  of  the  lateral  forces  is  about  an 
order  of  magnitude  less  than  that  of  the  normal  force  during 
buckling  events,  the  relative  change  in  lateral  forces  during 
buckle  formation  is  similar  to  the  changes  in  compressive 
force.  In  fact,  a  lateral  force  of  up  to  10  mN  was  observed  dur¬ 
ing  the  rapid  onset  of  buckling  for  the  1200  pm  array.  Addi¬ 
tionally,  forces  in  the  x-  and  y-orientations  may  change 
from  in  phase  to  out  of  phase  relative  to  each  other,  consis¬ 
tent  with  a  change  in  physical  orientation  of  CNTs  within  col¬ 
lapsed  cells.  SEM  examination  of  the  orientation  of  CNTs 
within  the  buckling  folds  often  reveals  a  preferential  shift  in 
lateral  direction  that  is  inconsistent  with  the  orientation  of 
the  CNTs  outside  of  the  indentation  region  (see  inset  of 
Fig.  6a).  The  localized  CNT  orientation  between  buckles  may 
change  by  as  much  as  90°  from  buckle  to  buckle.  The  contri¬ 
bution  of  lateral  forces  has  not  yet  been  explored  relative  to 
CNT  buckling  mechanics.  However,  based  upon  the  magni¬ 
tude  of  these  forces,  the  lateral  contribution  to  the  CNT  buck¬ 
ling  and  folding  should  not  be  overlooked  and  may  provide 
additional  evidence  of  cellular  collapse  consistent  with 
open-cell  foam  behavior. 
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Fig.  8  -  Lateral  and  normal  forces  measured  from  (a)  650  and  (b)  1200  jim  CNT  array. 


3.3.  Compressive  elastic  modulus 

The  compressive  elastic  modulus  of  the  CNT  arrays  along 
the  axial  direction  was  evaluated  using  the  CSM  technique 
simultaneously  with  force-displacement  data.  The  modulus 
values  reported  consider  full  contact  between  the  indenter 
tip  and  the  CNT  arrays  and  do  not  account  for  the  variation 
in  volume  fraction  of  CNTs  or  the  potential  increase  in  con¬ 
tact  area  between  the  CNTs  as  the  indenter  tip  load  is  in¬ 
creased.  These  assumptions  are  not  strictly  accurate,  as 
the  arrays  are  clearly  not  continuous  materials  but  are, 
rather,  composed  of  individual  CNTs  having  some  occupa¬ 
tion  fraction  much  less  than  unity.  Additionally,  the  free  tips 
of  the  array  likely  conform  to  the  indenter  face  at  greater 
loads,  increasing  real  surface  area  contact.  Nevertheless, 
the  assumption  simplifies  the  computation  of  material  prop¬ 
erties  and  allows  for  straightforward  comparison  of  results, 
particularly  in  the  absence  of  knowing  the  evolution  of  real 
contact  area  as  a  function  of  applied  load.  The  sample  elas¬ 
tic  modulus  was  then  computed  using  the  following  rela¬ 
tions  for  a  flat  punch  indenter: 


F  _  S 

r  2  Va 

1  _(i-y?)  |  (i-y?) 

Er  Es  Ej 


(5) 

(6) 
(7) 


where  S  is  the  unloading  stiffness  evaluated  at  40%  unload¬ 
ing,  AP  and  Ah  are  the  changes  in  applied  load  and  displace¬ 
ment  from  maximum  load  to  40%  unloading,  respectively,  E 
is  elastic  compressive  modulus,  A  is  cross-sectional  area  of 
the  indenter,  y  is  Poisson’s  ratio,  and  the  subscripts  r,  s,  and 
i  refer  to  reduced,  sample,  and  indenter,  respectively.  The  re¬ 
duced  modulus  is  that  which  is  directly  calculated  by  unload¬ 
ing  stiffness  during  indentation.  Poisson’s  ratio  is  a  very  small 
quantity  for  CNT  arrays  and  assumed  to  be  zero  for  computa¬ 
tion,  as  is  a  common  assumption  [28,34]  and  is  consistent 
with  our  experimental  observations. 

Fig.  9a  and  b  illustrate  the  average  compressive  elastic  mod¬ 
ulus  for  each  array  as  a  function  of  both  strain  and  raw  inden¬ 
tation  depth,  respectively.  Test-to-test  variation  in  modulus 
values  for  a  given  array  was  low  and  similar  in  magnitude  to 


Fig.  9  -  Average  elastic  modulus  of  CNT  arrays  as  a  function  of  (a)  strain  and  (b)  raw  displacement. 
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the  noise  range  of  the  measurement  for  these  particular  mate¬ 
rials.  Near  the  surface  of  the  CNT  arrays,  the  elastic  modulus 
monotonically  increased  as  contact  was  established  between 
the  tip  and  asperities.  Because  of  the  inherent  nonuniformity 
of  the  CNT  array  height  at  the  free  surface  due  to  asperities, 
the  initial  area  in  contact  with  the  indenter  face  was  likely  low¬ 
er  than  that  experienced  upon  further  penetration  into  the  ar¬ 
ray.  Thus,  the  reported  modulus  is  likely  a  conservative 
representation  for  the  initial  contact  region.  The  modulus  con¬ 
tinuously  increased  with  strain  to  a  point  approximately  corre¬ 
sponding  with  the  yield  plateau  of  the  stress-strain  profile. 
The  modulus  then  slowly  increased  with  depth  until  the  den- 
sification  point,  wherein  the  modulus  rapidly  increased.  Each 
of  the  arrays  in  this  investigation  exhibited  a  compressive 
modulus  value  between  10  and  20  MPa  at  the  initiation  of  yield. 
Remarkably,  the  magnitude  and  evolution  of  elastic  compres¬ 
sive  modulus  as  a  function  of  strain  is  quite  similar  for  all  CNT 
array  heights.  Because  the  deformation  is  largely  plastic  in 
nature,  the  consistent  modulus  values  may  be  a  result  of  elas¬ 
tic  recovery  of  the  top-most  asperities,  which  appear  relatively 
unaffected  by  indentation,  or  from  elastic  recovery  of  CNTs  be¬ 
neath  the  buckling  region. 


4.  Conclusions 

The  axial  compressive  mechanical  behavior  of  CNT  arrays 
with  heights  spanning  more  than  an  order  of  magnitude 
was  examined  to  significant  strains  using  flat  punch  nanoin¬ 
dentation.  It  was  found  that  the  compressive  mechanical  re¬ 
sponse  of  the  CNT  arrays  behaved  similarly  to  that  of  open¬ 
cell  foam  of  various  densities.  All  arrays  exhibited  initializa¬ 
tion  of  yielding  behavior  within  a  strain  of  0.3-0.12,  depend¬ 
ing  upon  the  array  height.  Well  organized  collective 
buckling  was  reported  for  the  first  time  originating  from  the 
top  surface  of  CNT  arrays,  with  an  additional  isolated  buckle 
near  the  substrate  surface  observed  for  arrays  obtaining 
heights  of  650  and  1200  pm.  This  localized  buckling  was  ex¬ 
plained  in  terms  of  the  local  cellular  geometry  variation  as  a 
function  of  depth  and  compaction  of  cellular  geometries 
within  the  arrays.  Lateral  forces  were  measured  for  the  first 
time  during  CNT  buckle  formation  and  indicate  a  significant 
lateral  component  of  buckling,  providing  additional  informa¬ 
tion  regarding  the  cellular  collapse  mechanisms.  The  magni¬ 
tude  and  relative  evolution  of  the  elastic  compressive 
modulus  as  a  function  of  strain  were  found  to  be  nearly  inde¬ 
pendent  of  CNT  array  height.  These  results  provide  significant 
new  insights  into  CNT  array  mechanical  properties  and  foam¬ 
like  behavior  as  a  function  of  array  height. 
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